Phase transitions of H2 adsorbed on the surface of single carbon nanotubes 
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By means of Diffusion Monte Carlo calculations, we obtained the complete phase diagrams of 
H2 adsorbed on the outer surface of isolated armchair carbon nanotubes of radii ranging from 3.42 
to 10.85 A. We only considered density ranges corresponding to the filling of the first adsorption 
layer in these curved structures. In all cases, the zero-temperature ground state was found to be 
an incommensurate solid, except in the widest tube, in which the structure with lowest energy is 
an analogous of the y/S x y/3 phase found in planar substrates. Those incommensurate solids result 
form the arrangement of the hydrogen molecules in circumferences whose plane is perpendicular to 
the main axis of the carbon nanotube. For each tube, there is only one of such phases stable in 
the density range considered, except in the case of the (5,5) and (6,6) tubes, in which two of these 
incommensurate solids are separated by novel first order phase transitions. 
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I. INTRODUCTION 



A carbon nanotube is a cylindrical structure^ that can 
be thought as the result of the wrapping up of a sin- 
gle graphene sheet^''^ over itself. In the same way that , 
in the past , graphene sheets were always found stacked 
to form graphite, carbon nanotubes were always found 
forming nanotube bundles, structures in which the only 
surfaces available for adsorption were the exposed parts 
of the tubes in the fringes of the bundle'^"^ . The adsorbed 
phases in a planar structure such as graphite are quite 
different from the ones observed on a highly patterned 
substrate like the external surface of a nanotube bundle. 
This can be seen by comparing the results on a bundle of 
carbon nanotubes and graphite for "^Hc^'^ (bundles)^'^° 
(graphite), Ha"'^^ (bundles)""!^ (graphite), and Ne^^'" 
(bundles) and^^'^^ (graphite). Similar to the opportu- 
nities for adsorption that the achievement of a single 
graphene sheet offers, a recent experimental work shows 
the proper technique to suspend a single carbon nanotube 
and study the phases of different gases (Ar,Kr) adsorbed 
on its external surface^'^. In this latter work, the results 
indicated that the phase diagrams were similar to those 
of the same gases on graphite, but with phase transi- 
tions shifted to higher pressures. In the present work, we 
have studied the adsorption behavior of H2 on a series of 
armchair ((n,n)) carbon nanotubes with different diame- 
ters in order to see if a quantum gas would exhibit novel 
phases when adsorbed on curved surfaces. The carbon 
nanotubes considered, together with their indexes and 
radii, are given in Table I. The most noticeable results 
driven from our calculations are: i), the existence of a 
solid-solid phase transition for (5,5) and (6,6) nanotubes, 
and ii), the stability of the curved version of the VS x VS 
commensurate phase only for the widest nanotube stud- 
ied (16,16). 



II. METHOD 

Our calculations were performed with the diffusion 
Monte Carlo (DMC) method. This well-established tech- 
nique allows us to solve exactly the Schroedinger equa- 
tion within some statistical errors to obtain the ground 
state energy of a system of interacting bosons^-'^. This is 
exactly our case, since in its lowest energy state (para- 
H2) the hydrogen molecule behaves as a boson. To apply 
the DMC algorithm, we need the potentials to describe 
the H2-H2 interaction (we chose the Silvera and Goldman 
expression^^, a standard model for calculations involving 
para-li2), and the C-II2 one. For the latter, we used 
a Lennard-Jones potential with parameters taken from 
Ref. 23, as used previously to describe the phase dia- 
gram of II2 adsorbed on graphene^'^. All the individual 
C-H2 pairs were considered, i.e., we took into consider- 
ation the corrugation effects due to the real structure of 
the nanotube. 

The last ingredient needed to perform a calculation in 
a DMC scheme is a trial wave function. It regulates the 
Monte Carlo importance sampling, and can be considered 
as a variational approximation to the exact description 
of the system. In this work, we used a trial wave function 
formed by the product of two terms. The first one is 
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where ri, r2, ..., r^ are the coordinates of the H2 
molecules and rj the position of the carbon atoms in 
the nanotube. The first term in Eq. (1) is a two-body 
Jastrow function depending on the II2 intermolecular dis- 
tances rij, while the second one has the same meaning 
but for each rij (C-H2). Finally, the third term is a prod- 



TABLE I. Armchair carbon nanotubes considered in this 
work, together with their tube radii (r^) and the most proba- 
ble distance for an adsorbed hydrogen molecule to the center 
of the cylinder (ro). In the last column, the adsorption energy 
of a single H2 molecule (eo) is reported. 

"Tube n (A) r-u (A) ep (K) 



(5,5) 

(6,6) 

(7,7) 

(8,8) 

(10,10) 

(12,12) 

(14,14) 

(16,16) 



3.42 
4.10 
4.75 
5.45 
6.80 
8.14 
9.49 
10.85 



6.36 
7.05 
7.70 
8.39 
9.75 
11.10 
12.46 
13.80 



-320.6 ± 0.1 
-329.0 ± 0.1 
-335.5 ± 0.1 
-341.0 ± 0.1 
-349.2 ± 0.2 
-353.7 ± 0.2 
-356.5 ± 0.1 
-357.2 ± 0.2 



uct of one-body Gaussians with n = \/xf+yf 

pcnd on the distance of each H2 to the center of the tube 

(ro)- 

Equation 1 is appropriate for describing liquid phases, 
but if the system is a solid, we need to multiply the trial 
wave function above (1) by 

JJexp {-C [{Xi - Xsitc)^ + iVi - 2/sitc)^ + {Zi - Zsito)^] } , 

i 

(2) 
whose purpose is to limit the location of the hydrogen 
molecules to regions close to the Xgite , 2/sitc , 2^sito regu- 
larly distributed coordinates of a solid arrangement. The 
variational parameters appearing in the whole trial wave 
function (&H2-H2: ^C-H2, ^ and c) were fixed to be the 
same values than the ones for H2 adsorbed on graphene^^, 
after having checked that there were no appreciable dif- 
ferences when variationally optimized for the (5,5) tube, 
i.e. 6H2-H2 == 3-195 A, 6c-H2 = 2.3 land a = 3.06 A^^. 
Parameter c was varied with density in the same way 
than for graphene. The only remaining parameter, ro, 
was optimized independently for each tube; the results 
are listed in Table I. 



III. RESULTS 

Table I contains also the binding energies (cq) of a sin- 
gle H2 molecule on top of the series of nanotubes consid- 
ered in this work. We observe, quite reasonably, that this 
energy increases with the radius of the cylinder. How- 
ever, even the value for the (16,16) tube is quite far from 
the result on a flat graphene sheet (-431.79 ± 0.06 K^'^), 
due to the surface curvature, that distorts the C-H2 dis- 
tances with respect to those of a flat graphene sheet. 

Our aim in the present work is to describe the pos- 
sible phases of H2 adsorbed on the surface of different 
(n,n) nanotubes. To do so, we considered the curved 
counterparts of all the commensurate solid phases found 
experimentally for most of the quantum gases ("^He, H2 
and D2) on graphite^'^^'^^'^®, which included the curved 
version of a a/3 x ^/3 phase, perfectly possible in these 
(n,n) tube substrates^^. For the incommensurate phases. 
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FIG. 1. Projection on a flat surface of an incommensurate 
structure corresponding to an areal density of 0.076 A~^ to 
be wrapped around a (5,5) tube. It can be described as an set 
of five intertwined helices whose pitch (indicated by an arrow) 
is 17 A. Every turn of the helix contains 10 H2 molecules. 



we tried structures similar to the triangular phases found 
in graphene, but wrapped up around to form hydrogen 
cylinders of radius ro (see Table I) . One of these struc- 
tures is shown in Fig. 1, which displays the projections 
on a plane of the H2 site locations for an incommensurate 
solid wrapped around a (5,5) tube. Here, the r coordinate 
represents the hydrogen positions on a circumference of 
radius ro = 6.36 A, while the z axis is chosen parallel to 
the main axis of the tube. One can see that the solid is 
built by locating five H2 molecules in circumferences on 
planes perpendicular to z, and rotating the molecules in 
neighboring circumferences half the distance between H2 
molecules. We defined a phase by the number of hydro- 
gen molecules in one of such circumferences, and varied 
their density by changing the distance between them. Al- 
ternatively, this structure can be thought as the result of 
having five helices of pitch 17 Awrapped around the tube, 
each one of them including 10 molecules per turn of the 
helix. 

DMC results for the equations of state corresponding 
to different phases of H2 adsorbed on a (5,5) tube are 
shown in Fig. 2. Open squares indicate a liquid phase 
(obtained by considering as a trial function only Eq. (1)), 
while the circles represent solid incommensurate phases 
with five (full circles) and six (open circles) H2 molecules 
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FIG. 2. Energies per H2 molecule for different adsorbed 
phases on a (5,5) tube. We display here the liquid (open 
squares), the different incommensurate solids: five- in-a- row 
(solid circles), six-in-a row (open circles), four-in-a-row (open 
triangles) and seven-in-a-row (full triangles). 



FIG. 3. Equation of state for the incommensurate solids ad- 
sorbed in carbon nanotubes of increasing radii: (5,5), (open 
squares, five-in-a-row solid; full squares, six-in-a-row arrange- 
ment); (6,6), (open circles, six-in-a-row structure; full circles, 
seven-in-a-row solid); (7,7), (open triangles); (8,8), (full tri- 
angles); (10,10), (inverted triangles). See further explanation 
in the text. 



per row. We display with triangles results for a four 
molecules per row (open), and seven molecules per row 
(full) arrangements. For all the areal densities consid- 
ered (up to 0.0937 A~^, the experimental density for a 
second layer promotion in planar graphite^®), these last 
two phases are unstable (of higher energy) with respect 
to the first two. All the alternative commensurate struc- 
tures are also unstable, as can be seen in Table II and 
III. In that table, we also display their areal densities, 
different for different tubes, even though the pattern of 
the adsorbed hydrogen molecules is the same. This is 
due to the fact that the areal densities depend on ro, 
while these phases are registered with respect to struc- 
tures whose dimensions depend on rt . Other solid incom- 
mensurate structures, such as helices with different num- 
ber of molecules per turn but with the same pitch, i.e., 
more or less H2 molecules on top of the dashed lines in 
Fig. 1, have energies greater than the phases represented 
by the circles in Fig. 2. For instance, a structure with 
eleven molecules per turn, instead of the ten displayed 
in Fig. 1, has a minimum energy of -336.0 ± 0.1 K, for 
a density of 0.062 ± 0.003 A-^. Therefore, the ground 
state of H2 adsorbed in the outer surface of a (5,5) car- 
bon nanotube is an incommensurate solid with five Il2's 
per row. Upon a density increase, there is a first-order 
solid-solid phase transition to another incommensurate 
solid, similar to the first one but with six atoms per row. 
The equilibrium densities for both structures at the tran- 
sition are obtained from a double tangent Maxwell con- 
struction: 0.0685 ± 0.0001 A~^ (five molecules per row, 
energy per hydrogen molecule -345.3 ± 0.1 K) and 0.0795 
± 0.0001 A~^ (six intertwined helices, with energy per 



hydrogen molecule -334.7 ± 0.1 K). 

The stable phases of H2 on other (n,n) tubes, up to n = 
14, are similar to the (5,5) one. In all cases, the ground 
states are the same type of incommensurate solids al- 
ready described, the only difference being the number 
of intertwined helices (or molecules in the same circum- 
ference) forming the structure. This number was found 
to be always equal to the index n of the nanotube. In 
Fig. 3, we show the equation of state for some of such 
solids, from the already displayed (5,5) case (open cir- 
cles; five molecules in a row), to the (10,10) one (inverted 
open triangles; ten molecules in a row), going through the 
(6,6) (open circles), (7,7) (open triangles) and (8,8) (full 
triangles) tubes. In all cases the dashed lines are mere 
guides-to-the-eye. Between the two narrowest cylinders 
and the rest there is an important difference, though: the 
(6,6) tube exhibits a first order phase transition between 
two incommensurate solids with six and seven molecules 
per row while, while for the other tubes, there is only 
one stable incommensurate structure in the areal den- 
sity range corresponding to an adsorbed monolayer The 
zero-pressure densities and energies corresponding to the 
ground state for all the tubes considered are given in 
Table II. There, and in Table III, we can see that those 
incommensurate solids are the truly ground states for the 
systems under consideration at zero pressure, since their 
energies are lower than the corresponding to any of the 
commensurate structures on the same tubes. The only 
exception is the (16,16) nanotube, in which the ground 
state is the same v3 x vS structure than in a flat sur- 



TABLE II. Energies per particle (e;,) and equilibrium densities (p) for different phases proposed in the 
gases on graphite, when adsorbed in tubes of different radii. The hquid and incommensurate helical 
corresponding to the minimum energies obtained by means of third-degree polynomial fits to curves 
Fig. 1. The rest correspond to exact densities. Error bars within parenthesis represent the uncertainty 



literature^ for quantum 
densities are the values 
of the type displayed in 
of the last figure shown. 
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liquid 




2/5 




3/7 




73 X \/3 




helical incommensurate 
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P(A-^) 


e, (K) 


P (A-^) 


eb (K) 


P (A-^) 


eb (K) 


P (A-^) 


et, (K) 


P (A-^) 


et (K) 


(5,5) 


0.0574(4) 


-340.3(1) 


0.0407 


-333.7(2) 


0.0436 


-325.4(2) 


0.0339 


-329.03(4) 


0.0621(3) 


-349.0(2) 


(6,6) 


0.0586(3) 


-349.70(9) 


0.0441 


-344.4(1) 


0.0472 


-334.3(2) 


0.0367 


-339.52(6) 


0.0673(1) 


-360.34(4) 


(7,7) 


0.0615(6) 


-356.2(2) 


0.0471 


-352.28(7) 


0.0504 


-339.7(2) 


0.0392 


-347.66(3) 


0.0716(1) 


-367.89(5) 


(8,8) 


0.0611(7) 


-360.3(2) 


0.0494 


-358.9(1) 


0.0529 


-343.4(2) 


0.0411 


-354.44(3) 


0.0742(1) 


-372.43(3) 


(10,10) 


0.0578(3) 


-365.95(6) 


0.0531 


-367.11(8) 


0.0569 


-347.5(2) 


0.0443 


-364.12(5) 


0.0758(1) 


-375.6(1) 


(12,12) 


0.0565(4) 


-369.1(1) 


0.0559 


-371.8(1) 


0.0600 


-346.3(3) 


0.0466 


-369.93(5) 


0.0703(2) 


-375.96(5) 


(14,14) 


0.0530(4) 


-370.4(1) 


0.0581 


-373.5(1) 


0.0623 


-342.5(3) 


0.0485 


-373.10(8) 


0.0681(2) 


-375.3(1) 


(16,16) 


0.0511(8) 


-367.57(9) 


0.0600 


-371.17(5) 


0.0643 


-333.8(3) 


0.0500 


-371.7(1) 


0.0660(1) 


-371.4(1) 



TABLE III. Same than in Table II, but for two other commen- 
surate structures found experimentally for D2 on graphite. 
The dimensions of their unit cells make them only possible 
for the tubes shown. In all cases, they are unstable with re- 
spect to the incommensurate structures (n< 14) or to the 
\/3 X ^/3 one (n= 16). 



T/ie^"-^" 


Tube 




P(A-^) 


et, (K) 


(8,8) 
(12,12) 
(16,16) 




0.0540 
0.0612 
0.0656 


-359.4 ± 0.1 
-371.6 ± 0.1 
-369.9 ± 0.1 


31/75 (5 phase in H2 on graphite^"') 


Tube 




P 


Eb 


(5,5) 
(10,10) 




0.0421 
0.0549 


-334.8 ± 0.1 
-367.0 ± 0.1 



face. The stability limits for the two solids in the (6,6) 
tube are 0.0755 ± 0.0001 A^^ (et = -338.2 ± 0.1 K) and 
0.0855 ± 0.0001 A-2 (cb = -354.1 ± 0.1 K), also ob- 
tained by a Maxwell construction. The 1 x \/3 structure 
defined in Ref. 24 for (n,0) nanotubes was found also to 
be unstable for all the tubes considered. For instance, 
for the (16,16) tube, the density of this phase is 0.0750 



A ^, with a binding energy of -303.9 ± 0.6 K. 



IV. CONCLUDING REMARKS 

Summarizing, we studied all the possible H2 phases 
adsorbed on the surface of armchair carbon nanotubes 
ranging from (5,5) to (16,16). In all cases, but the last 
one, we have found that the stable phases are incom- 
mensurate solids formed by molecules adsorbed on cir- 
cumferences perpendicular to the main tube axis. Where 
the curvature of the surface is more relevant, i.e., in the 
narrowest tubes [(5, 5), (6, 6)], our results show the exis- 
tence of a solid-solid zero-temperature phase transition 
between two incommensurate structures. The first com- 
mensurate solid phase, curved version of the well-known 
v3 X v3 phase, appears only when the radius of the nan- 
otube is big enough: a (16,16) tube for H2. 
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